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Andreev scattering and the Josephson current through a one-dimensional interacting electron 
liquid sandwiched between two superconductors are re-examined. We first present some apparently 
new results on the non-interacting case by studying an exactly solvable tight-binding model rather 
than the usual continuum model. We show that perfect Andreev scattering (i.e. zero normal scat- 
tering) at the Fermi energy can only be achieved by fine-tuning junction parameters, a fine-tuning 
which is possible even with bandwidth mismatch between superconductor and normal metal. We also 
obtain exact results for the Josephson current, which is generally a smooth function of the supercon- 
ducting phase difference except when the junction parameters are adjusted to give perfect Andreev 
scattering, in which case it becomes a sawtooth function. We then observe that, even when interac- 
tions are included, all low energy properties of a junction {E « A, the superconducting gap) can 
be obtained by "integrating out" the superconducting electrons to obtain an effective Hamiltonian 
describing the metallic electrons only with a boundary pairing interaction. This boundary model 
provides a suitable starting point for bosonization/renormalization group/boundary conformal field 
theory analysis. We argue that total normal reflection and total Andreev refiection correspond to 
two fixed points of the boundary renormalization group. For repulsive bulk interactions the An- 
dreev fixed point is unstable and the normal one stable. However, the reverse is true for attractive 
interactions. This implies that a generic junction Hamiltonian (without fine-tuned junction param- 
eters) will renormalize to the normal fixed point for repulsive interactions but to the Andreev one 
for attractive interactions. An exact mapping of our tight-binding model to the Hubbard model 
with a transverse magnetic field is used to help understand this behavior. We calculate the critical 
exponents, which are different at these two different fixed points. 



I. INTRODUCTION 

One of the fascinating consequences of superconductivity is the phenomenon of Andreev scattering at a normal 
metal/superconductor (NS) junction. This corresponds to an incoming electron from the normal side being reflected 
back as a hole, thereby producing an additional Cooper pair in the superconductor condensate. Both normal and 
Andreev reflections are expected to occur in a realistic NS junction, in addition to quasiparticle transmission into the 
superconductor. If the gap is large enough, the latter's propagation is suppressed and squared reflection amplitudes 
add up to one through probability current conservation. Building up on Andreev reflection, one arrives at the related 
phenomenon of the Josephson current, in which the normal region of an SNS junctions carries a supercurrent driven 
by the gap phase difference between left and right superconductors. 

While the original work in this fieldElu treated the normal metal within the Fermi liquid framework {i.e. essentially 
ignored interactions), the effect of interactiiQiis for the case of a one-dimensional metal between two superconductors 
has been treated recently by several groupsfl B using bosonization and renormalization group methods. The methods 
and conclusions of these works are closely related to previous work on tunnelling through a single impurity in a 
quantum wireB 

We have chosen to re-examine this subject in both non-interacting and interacting cases because we feel that 
previous treatments have missed some interesting physics. In particular, most of the standard work on the non- 
interacting case has essentially ignored band structure effects, using a free electron model with a pairing potentijal 
that varies abruptly accross the junction, together with a scattering potential at the interface. A more recent papeit£l 
considers the case where the Fermi velocity is different on the S and N side. The conclusion of this work is that, at the 
Fermi energy, there is perfect Andreev reflection (and therefore normal reflection) when the scattering potential and 
velocity mismatch are absent. Both of these effects serve to increase the normal scattering amplitude in an additive 
way. We arrive at qualitatively different conclusions by explicitly including band structure in the form of an exactly 
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solveable one-dimensional tight-binding model with Hamiltonian: 

H-HN = J2 [i-tj,3+l^Ja^j+la + AjVltV'ji + h.C.) + {Vj - • (1-1) 

j 

The chemical potential, fi, is assumed to lie within the band on the normal side \fi\ < 2t and h.c. stands for Hcrmitcan 
conjugate. Here the interface is chosen to lie between sites and 1 with the superconductor on the negative x-axis 
and the normal metal on the positive x axis so that: 




Vi=VSiU A,= (1.2) 



Aj represents the pairing interaction which exists on the superconducting side (j < 0) only. For simplicity we consider 
both normal and superconducting sides to be one-dimensional, but see below. The notion of a "perfect junction" 
becomes less clear in such a model. In the particle-hole symmetric case, jj. = V = 0, we find that there is always 
some normal scattering at the Fermi energy unless the interface tunnelling parameter, t", is fine-tuned to a particular 
value. In the limit |A| «t,t', this particular value becomes 

t" = Vtt'. (1.3) 

For general values of the chemical potential we find that both t" and the normal scattering intensity, V, must be fine- 
tuned in order to achieve perfect Andreev reflection. Wc emphasise that these conclusions seem to be different than 
the previous ones obtained without explicit consideration of band effects. For example, in the particle-hole symmetric 
case it is possible to get perfect Andreev reflection even with Fermi velocity mismatch {t ^ t') provided that t" is 
adjusted to the right value. It is worth emphasizing that the case of infinite interface tunnelling: t" » t,t' does not 
correspond to perfect Andreev reflection as one might naively suppose, but instead to zero Andreev reflection. The 
physical reason is that, in this limit, two electrons get trapped at the interface on sites and 1, effectively decoupling 
all sites with j < from all sites with j > 1. In this limit the normal side does not "feel" the pairing and hence 
exhibits no Andreev reflection. 

We arc not aware of any previous explicit calculation of the Josephson current for such an interface model which we 
perform here. We consider two, possibly different, interfaces separated by a distance of I lattice sites with the pairing 
potentials having a phase difference x- The (zero temperature) Josephson current is defined from the derivative of 
the groundstate energy with respect to this phase difference: 

/(X) = 2e^Eo. (1.4) 
ax 

While the groundstate energy is obtained by summing over all states below the Fermi surface, we show explicitly 
that its derivative with respect to x only depends on quantities deflned at the Fermi surface, being insensitive to the 
details of the band structure. When the junction parameters are fine-tuned to give perfect Andreev reflection we find 
that the Josephson current is a sawtooth function of x- 

lix) - ^X, (mod 27r), (1.5) 

with steps of size 2evf/l occurring at x = (2n + 1)tt. For any other choice of junction parameters the Josephson 
current is a smooth function of x- For example, in the limit of a weak junction, t" « t,t' , we find: 

7(x)cx(i")'sinx. (1.6) 

While these formulas are well-known results, our method in fact allows us to calculate exactly the zero-temperature 
Josephson current in the general case of arbitrary amounts of normal versus Andreev refiection, independently at each 
boundary. The full crossover between the above two limits is thus described. 

One approach which we shall make much use of in this paper is that of an effective field theory, whereby the 
superconducting side is replaced by a particular boundary contribution on the normal side. There are several reasons 
why it is convenient to "integrate out" the electrons on the superconducting side of the junction in such a way as to 
obtain an effective Hamiltonian for the electrons on the normal side. This is a very natural thing to do considering 
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the fact that the superconducting electrons have a gap in their spectrum: if we consider physics at energy scales small 
compared to the gap we expect to obtain a simple effective action without any retarded interactions. (This may break 
down for non s-wave pairing where the gap vanishes in certain directions; we do not consider that case here.) The 
resulting effective Hamiltonian (for a single junction) is: 

H-^iN = -tJ2 [(V']aV'j+i<T + h.c.) - fii'j^i'ja] + [AbVI,t^u + '^■^] + VB^pL^pia- (1.7) 
The effective boundary pairing interaction, As and effective boundary scattering potential, Vb, depend on all the 



parameters of the superconductor and the junction, t', A, t" , V . Beginning from our interface model of Eq. (1.1) 
we determine explicitly the parameters As and Vb of the boundary model and check that low energy properties are 
faithfully reproduced. Of course, we again find with the boundary pairing model that perfect Andreev reflection only 
occurs if the boundary parameters are fine-tuned. In particular, for the particle-hole symmetric case, /i — Vb = 0, 
the condition is: 

\^B\=t. (1.8) 

Note that is is not Ab — > oo as one might naively suppose. 

One advantage of the boundary model is that it should arise from much more general, and more realistic, interface 
models. While the simple form of Eq. ( [l.l| ), quadratic in fermion operators, is the result of a mean field approximation 
to a more realistic model with pairing or electron-phonon interactions on the superconducting side, we expect that the 
effective Hamiltonian of Eq. (^]^) will still be valid at energies small compared to the gap when the interactions on 
the superconducting side are treated more accurately. Furthermore, it is more or less obvious that the same effective 
Hamiltonian arises when the one-dimensional normal metal is coupled to a three-dimensional superconductor. This 
is an important generalization since superconductivity is not believed to occur in a strictly one-dimensional system. 

More generally we wish to add interactions to our Hamiltonian on the normal side which we assume to be one- 
dimensional. A simple choice would be an onsite (Hubbard) interaction for all sites j > 1- 



H ^ H + Hint 

2 V 



H^nt = ^Y.^n, - 1)\ (1.9) 

where Uj is the total electron number operator at site j. We could also consider longer range density-density inter- 
actions. We will be interested in the case of both repulsive and attractive bulk interactions; the latter may arise in 
a low energy effective theory from phonon exchange. The negative U Hubbard model has a gap for spin excitations 
which can be eliminated by considering longer range interactions. We will discuss both cases with zero and non-zero 
spin gap. These interactions can be treated essentially exactly, at low energies, using bosonization, renormalization 
group and conformal field theory techniques. It is a major purpose of this paper to discuss how to generalize these 
techniques to the interface model. We a rgue that this is best done by integrating out the superconducting electrons to 



obtain the boundary model of Eq. (1.7) with the bulk interactions. Hint added. (More generally, we might also obtain 
additional interactions at the boundary. These can be treated in the same framework.) Indeed, this approach seems 
to be more or less forced upon us by the renormalization group philosophy of integrating out high energy modes to 
obtain an effective low energy Hamiltonian. Having performed the initial step of integrating out the gapped degrees 
of freedom on the superconducting side we may then proceed to analyse the effective Hamiltonian with bulk and 
boundary terms using general methods developed to deal with quantum impurity problemsO Of course, our results 
will only be valid at energies E << A. 

The boundary renormalization group approach leads to the conclusion that the boundary interactions will renor- 
malize to a fixed point corresponding to a conformally invariant boundary condition. It appears likely that there are 
only two such boundary conditions that occur in this problem, in the particle-hole symmetric case (ii = V = Vb = 0) 
corresponding to a free boundary condition (b.c), Ab ~ which preserves electron number and therefore has no 
Andreev reflection and to an "Andreev boundary condition" for which there is perfect Andreev reflection. Which 
of these boundary conditions is stable under renormalization group transformations depends on the sign of the bulk 
interactions, U . We find that the free b.c. is stable for repulsive bulk interactions but the Andreev b.c. is stable 
for attractive bulk interactions (f/ < 0). We calculate the various critical exponents associated with these critical 
points. It is important to realize that critical exponents are characteristic of a particular fixed point and are different 
at the free and Andreev fixed points. Thus, for instance with repulsive bulk interactions, if we started with bare 
interface parameters that put the interface close to the Andreev fixed point then the exponents characterizing the 
initial flow away from the Andreev fixed point at high temperature are different than the exponents characterizing 
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the flow towards the free fixed point at low temperature. However, it should be emphasized that for a generic choice 
of bare interface parameters the Hamiltonian will not be near the Andreev fixed point; this requires fine-tuning. The 
situation is similar in the non particle-hole symmetric case except we now get lines of fixed points with either or 
perfect Andreev reflection. I-. 

The first work that we are aware of on Andreev scattering in Tomonaga-Luttinger liquidsO attempted to apply 
boundary RG techniques without explicitly intergrating out the superconducting side and without taking into account 
the effect of the b.c.'s on the exponents. This led to incorrect predictions for the exponent governing the Josephson 
current. [|^ used a method closely related to ours but applied it to a different geometry: a closed normal ring in contact 
with superconductors at two points. TakaneLl corrected some of the earlier result in [ [sj for the exponent governing 
the Josephson current using methods essentially equivalent to ours but without explicitly invoking the concept of 
integrating out the superconducting side. We extend Takane's result by introducing the conceptually important and 
very useful notion of integrating out the superconducting electrons thus making clear the relationship between the 
interface problem, other quantum impurity problems and boundary conformal field theory. This facilitates a more 
general discussion of the universal critical behaviour. In particular we discuss the behaviour of the Josephson current 
in the vicinity of the Andreev fixed point, obtaining quite different results than those in [ ^, and discussing how 
the functional dependence of the current on the superconducting phase difference crosses over between sawtooth and 
smooth forms. 

Given the difficulty of achieving perfect Andreev scattering in the non-interacting case our conclusion is quite 
remarkable that, with attractive bulk interactions, a generic interface will renormalize to perfect Andreev scattering 
as r 0. It must be admitted that this conclusion is based on an unproven but widely made assumption about RG 
flows and fixed points in the boundary sine-Gordon model which arises here after bosonization. In order to make some 
of our rather unintuitive results seem more plausible we discuss an exact mapping of our boundary pairing model 
(in the particle-hole symmetric case) into a Hubbard model with a bulk magnetic field and a transverse boundary 
magnetic field. In the transformed model the Andreev boundary condition corresponds to one in which the boundary 
electron has a frozen transverse spin polarization. 

In the next section we give the solution of the interface and boundary models and discuss their equivalence, in 
the case of zero bulk interactions. Most of the details of the equivalence are relegated to an appendix. In Sec. Ill 
we consider an S1NS2 system and calculate the Josephson current. In Sec. IV we include bulk interactions in the 
boundary model and determine phase diagrams and critical exponents. In Sec. V we discuss the exact mapping onto 
the Hubbard model with bulk and boundary magnetic fields. 



II. THE LATTICE INTERFACE AND BOUNDARY MODELS 



In this section, we consider various lattice models of normal me tal-superconductor contacts in the case of zero bulk 
interactions. The generic Hamiltonian for all of these will be (p~l|), with various geometries, sets of hopping strenths, 
and potentials specified along the way. The calculation procedure is similar to the usual one used in dealing with 
continuum models: the wavefunctions are found for each sector, and the matching conditions at the contacts yield 
consistency equations from which the various reflection/transmission coefficients are obtained. 

We start by performing the Bogoliubov-de Gennes transformation 

^J"'^ = ~ ^^ajli-a], (2-1) 

a 

where the quasiparticle operators satisfy {Tqct, = ^aa'Saa' and a is a (real) quasimomentum index. We obtainS, 

by requiring the Hamiltonian to be diagonal {H = Eq + J2a^°''^a(7laa) the following lattice Bogoliubov-de Gennes 
equations: 

eaVaj = + + + 1 - {Vj ~ n)Vaj + A*Uaj. (2.2) 

The solutions of these equations in two particular geometries are presented below. The emphasis is put on the 
calculation of the Andreev reflection coefficient and on the low-energy properties of the system. 



A. The lattice interface model 



As a simple toy model for a NS interface, we consider the geometry depicted in figure |l|. This tight-binding model 
is one of free electrons in which a nonvanishing superconducting order parameter has been induced by an unspecified 
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mechanism on the left-hand sites of the lattice only. Although it is a straightforward lattice version of the ubiquitous 
continuum one-dimensional models used in [ ^ and numerous subsequent work, let us underline that our model 
includes some additional features, namely different bandwidths t, t' on the normal and superconducting sides, and an 
arbitrary coupling t" together with a local scattering potential V at the interface. In other words, our lattice interface 



model is defined by the BdG equations ( ]2.2D with the choice of parameters in (1.2). 

A A A V 
••• • t'»t"« t • t • ... 

-2-1 12 3 

FIG. 1. The interface model. Dots represent lattice sites, labeled by integers. On the left is the superconductor, with gap A 
and bandwidth t' . It is coupled to a normal metal with bandwidth t through a junction of strength t" , with a local potential 
V to tune the normal scattering rate at the contact. 

The calculation of the normal and Andreev reflection coefficients in the lattice interface model is nothing but 
an exercise in elementary quantum mechanics. Namely, it is performed by choosing an appropriate ansatz for the 
wavef unctions, after which solving for the matching conditions of the eigenstates at the interface yields all the desired 
information. 

Let us implement this procedure by taking a simple travelling wave solution to the Bogoliubov-de Gennes equations 



(2.3) 



In the above equation, the quasimomcnta j3 and 5 will turn out to be complex for the region of parameters we will be 
concentrating on, namely at energies below the superconducting gap. This means that the wavefunction amplitudes 
Uaj and Vaj are exponentially decaying on the left-hand side. 



The rest of the procedure is then a simple matter of substituting (2.3) in (2.2) and carrying out the necessary 
algebra. The energy is 

Cq = — 2t cosa — (2.4) 
while the other parameters in (|2.3|) are given by 



C0S7 = cosa -|- /i/t, 2t' cos [3 ^ \J (2i cosa + /^)^ — |Ap — ^,, 

2t' cos 5 = ^J{2tcosa + nY - |A|2 -f- fx. (2.5) 

Although all reflection and transmission coefficients can be calculated explicitly, a simpler form for these expressions 
is obtained if we are interested primarily in the scattering of particles whose energy is very small compared to the 
superconducting energy gap. The e — > limit gives 

cos a = — cos 7 = ^ H J sin a = sin 7 = — /.i^ /At-^, 

cos f3 = -{cosSy* = -i\A\/2t' - fi/2t', sin/3 = (sin(5)* = yr~^i|A|TAi)V4t^, (2.6) 

and allows us to write the following expression for the Andreev reflection coefficient: 

^2 y/l - ;^2/4t2 [-1^ +smP + sin (5 



RaU^o = ie 



tt" 



(2.7) 



In the above, x is the phase of the order parameter, A = | A|e'^. Note that the reflection coefficients moreover obey the 
sum rule \Ra\'^ + I-Rat P = 1 at half-filling, in view of the conservation of the current tjj+i [u*^jj^iUaj — v*^j+iVaj — h.c.]. 

This low-energy Andreev reflection coefficient, as a function of the bandwidths t, t' and of the strength of the 
pairing A, as well as of the tunneling strength t" and local scattering potential V , obeys some simple but interesting 
properties. Imagining for example a generic experimental setup in which variations in i" and V can be implemented, 
we can ask how RA\e=o behaves. A simple variation yields that the optimal amplitude is achieved for 
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2 2t 

(Here c.c. denotes complex conjugate.) Subsequently tuning t" yields a maximum at 



^ max ~ I . I , . „ ■ (2-9) 



+ (sin/3 + c.c.) 
In the particle-hole symmetric case, this simplifies to 

t"'^ = tt' 

max '''' 



1 



lAI 



(2.10) 



4t'2 2t' 

Tuning the parameters in such a way, we find perfect Andreev reflection, i.e. 

i?A |e=o, ^^e-^'' (2-11) 
In terms of left- and right-movers in the continuum limit of ( |2.3| ), defined in Sec. IV, this corresponds to 

vE'flt(0) = -ze'^*[^(0), vl/^i(0) = ze'XvE'[^(0), (2.12) 

which are simply the perfect Andreev conditions usually imposed. 

Thus, pure Andreev boundary conditions are obtainable in the lattice interface model by simply tuning two param- 
eters in the general non particle-hole symmetric case. For the particle-hole symmetric case, it is enough to tune one 
parameter. 

B. The lattice boundary model 

At energies below the superconducting gap, all wavefunctions incident on the NS interface from the normal side are 
eventually reflected back through either normal or Andreev reflection processes. The vanishing of the transmission 
coefficients opens the door to the formulation of a different approach than that adopted in the interface model, namely 
one in which we consider a system with a boundary obtained by "integrating out" the superconducting side, yielding 
an effective pairing potential localized at the original contact. In all generality, we will also consider a potential 
present at the boundary. 

What we could call the lattice boundary model can be directly defined through the BdG equations ( |2.2D , but our 
system (illustrated in figure |^) is now taken to live on the axis of positive integers, with a boundary at j = having 
on-site pairing As and potential Vb on the first site. The hopping strength is taken to be t between all sites. For the 
moment, we will perform all the relevant calculations from scratch, deferring the explicit connection with the previous 
interface model until later. 



FIG. 2. The boundary model. The superconductor has been replaced by effective boundary potentials As and Vb- 

In the same way as for the interface model, we can calculate the normal and Andreev reflection coefficients by 
simple quantum mechanics. Of primary interest is the Andreev one, which is readily obtained by using the ansatz 

u^j- = e-"J' + RNe'"^, Vaj = (-l)^'i?Ae*^^' (2.13) 

in the lattice BdG equations ( ^.2| ) in the geometry just described. We find again — -~2tco8a — ji with cos 7 — 
cos a -f /i/t, together with 
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A^t sina /o t a\ 

^ ^ |AbP + (te— + yB)(-te-7 + Vb) ^ ' 



which, at zero energy, becomes 





FIG. 4. |_Ra| plotted against lAsl/t for Vs/t = 0, 1,2, and 4. 



The modulus of this coefficient is plotted in figure I at half-filling, as a function of the boundary pairing As and 
boundary scattering potential Vb- The maximal amplitude has unit modulus, and such a maximum occurs for any 
choice of filling in view of a similar optimization property to the one in the interface model. Namely, tuning Vb and 
Ab yields a maximum amplitude at 



and this once again produces the perfect Andreev condition 



(2.16) 



(2.17) 



Figure (|j) shown the influence of a progressively stronger boundary scattering potential Vb on the Andreev reflection 
coefficient at half-filling, plotted against the boundary pairing strength. 

The equivalence of formulas (2.7) and (2.15) for the effective low-energy Andreev reflection coefficients comes from 
the fact that, as mentioned before, the lattice boundary model is obtainable by integrating out the gapped side of the 
interface model. This procedure is outlined in Appendix A, in which the explicit relationship between interface and 
boundary pairings and potentials is derived. At low energies, this correspondence reads 



2V 



t' 



sin j3 - 



Vb = V + 



21/ 



i sin P + i sin 5 



(2.18) 
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In the particle-hole symmetric case, this means that 



l + |A|2/4t'^- |A|/2t' 



(2.19) 



The interesting aspect of this formula comes from its behaviour in various limits. Contrary to simple intuition, a 
large bulk pairing |A| >> t' on the superconducting side induces a small boundary pairing As, according to the limit 



(2.20) 



" |A|»t' |A|' 

whereas a small |A| << t' yields a boundary pairing whose amplitude depends on the hopping parameters exclusively: 



A, 



^e-"^ — 

|A|<t' t' 



(2.21) 



It is important to note that these formulas hold for e ^ |A|, so one should not be surprised that a small bulk pairing 
still produces a significant boundary pairing, with a finite amount of Andreev reflection. The limits of zero energy 
and zero pairing do not commute. 

The above formulas allow one to move freely between interface and boundary formulations of the NS problem, as 
long as the low-energy sector of the theory is considered. They will be used in the next section, which is devoted to 
the calculation of the Josephson current in a SNS junction. 



III. JOSEPHSON CURRENT 



From the considerations of the earlier sections, we see that the problem of an S1NS2 superconducting junction can 
be investigated within a double boundary framework, provided we are interested only in energies much smaller than 
the gaps on either side. If we imagine integrating out both the left and right superconductors, we obtain a lattice 
model with two boundary pairings, which we dub the boundary junction model. Namely, we take this to be the model 
of free electrons in the bulk, with pairings A^ and A^ at the right and left ends. It is important to realize that 
A_R,L are boundary pairings, whose influence on the reflection coefficients has been explained in detail in the previous 
section. One should be careful not to confuse t hem with bulk pairings, which have an altogether different effect (the 
two are related, of course, by the relationship ( 2.18| )). We thus define the system on sites 1, I — 1 and take 



Ai = Al, Ai_, = Ane'^, (3.1) 

with Aji i^ G R (that is, we have put all the superconducting phase difference on the right pairing). 
Solving the lattice BdG equations by using the ansatz 

Uaj — Aa sinaj -I- cosaj, Vaj = (—1)"' (Ca sinaj -I- Da cos aj) (3-2) 

with e = —2tcosa yields after a certain amount of algebra the condition for the allowed quasiparticle momenta a (we 
have chosen / to be odd). We find (for convenience, we have set i = 1 in what follows) 

= sin^ al + 2ArAl cos x sin^ « - {Ar + A^) sin^ a{l - I) + A%Al sin^ a{l ~ 2). (3.3) 

Using a generalization of the approach used in [ ^ to treat the problem of free electrons on a tight-binding chain 
with a boundary scattering potential, we can conveniently find the closed form solution. First of all, let us write the 
allowed momenta in terms of energy-dependent phase shifts d± as 

an± = — + — , (3.4) 



where ± refers to the two independent sets of Andreev levels, labeled by the integer n. Substituting this into (3.3) 
yields 

/i cos 2(5 sin 2(5 = /a /i= 1 - (A^ + A|) cos2a -I- A^A^ cos4a, 

/2= -{Al + Al) sin2a + A^Al sin 4a, 

/3= (1- A2j)(1- A2)+4AHAisin2acosx. (3.5) 
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Noting that /| < + /| throughout the parameter space allows us to write this as 

cos{2S — g) — cosh, (3.6) 

where the functions g and h are defined by 



cos.g = ^= cos h = (3.7) 



Studying carefully the various functions along paths in the parameter space yields a consistent choice of branches for 
the inverse trigonometric functions. This leads to the final answer for the phase shifts, which we write as 



where 



ig±h)/2 (3.8) 



-27r — g a < a , 

g a < a < TT — a g = sgn{a ~ 7r/2) arccos — 

2TT-gTT-a<a V /i + /; 



2 ' 



a = ^3? {arccos ^(A^^ + A^^)} , h — arccos ^ — - (3.9) 

V /r + 

and all arccos functions have their image in [0, tt]. The phase shifts are analytic except at the bottom and top of the 
band at the perfect Andreev points, where g suffers a branch jump. 

Knowing the phase shifts allows us to compute all the energy levels, and understand their behaviour as a function 
of the boundary pairings and the superconducting phase difference x- One recovers the usual picture wherein the 
levels are separated by finite gaps on an energy versus phase diagram, except at the perfect Andreev points, where 
the gaps vanish. 

The reason why it is so convenient to solve for the allowed quasiparticle momenta in terms of these phase shifts, 
is that this procedure allows us to write the ground-state energy straightforwardly in a l/l expansion. Again, the 
derivation is very similar to the one in [ and we refer the reader there for the missing details. 

The ground-state energy can be written as the sum of the individual energies of the occupied Andreev levels. 
Namely, 

N 

^o = EE^K±]- (3.10) 

n=l ± 

In the limit of large I, we can use a Euler-MacLaurin formula to transform this sum into an integral. Subsequently 
expanding to order l/l, we get 



Eo^2l — e(fc) + - / de{d+{e) + <5_(e)) , , 

Jo TT Jeo l 



1 f6+{kF)Y 1 fS^(kF)Y 1 



2 \ TT J 2 \ TT J 12 



(3.11) 



where kF = 7r(A^ -I- \/2)/1,€f = e{kF),£o = e(0) = —2,vf = e'{kF)- The crucial thing to notice here is that the l/l 
terms are functions of data exclusively at the Fermi surface. While the term is well-known to correspond to the 
finite-size contribution from open boundary conditions for a conformal field theory with central charge c — 2 like the 
present one (each spinful chiral fermion carries a unit conformal charge) , the other terms depending on the squares of 
the phase shifts at fcjr give the change in Eq coming from the effect of the boundary pairings. The 0(1^) term, given 
by Fumi's theorem, depends however on S across the whole filled part of the band. 

This expression for the ground-state energy allows us to write the Josephson current at zero temperature in the 
limit of large I, in the presence of arbitrary boundary pairings, i.e. with an arbitrary amount of normal versus Andreev 
reflection on either edge. Namely, the Josephson current is given by 

lix) = 2e^Eo. (3.12) 
ax 

Upon calculating this derivative, one easily sees that Fumi's theorem 0{l'^) term does not contribute to I{x), since 
the sum of the phase shifts is independent of x for any energy (in ( |3.8[ ), only h depends on x)- Thus, the Josephson 
current is controlled exclusively by parameters at the Fermi surface, and is given by the general expression 
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I{X) = -^^R^L 





g(cosx) 




sin X arccos 


[(g2(i)+A2t,|,)i/2j 





4AflAi sin2 |5(cos2f) + A2 



(3.13) 



where we have defined 

~g{y) = (1 - A2,)(1 - Ai) + A^A^^y, 



A= (Afl- AL)(l + Ai?Ai). 



This function reproduces the well-known behaviours in the limiting cases of perfect or very weak Andreev reflection: 
the fine-tuning for perfect Andreev reflection on both sides of the junction corresponds to setting A^j = A^ = 1, 
which, when substituted in (3.13), yields 



(Ishii's sawtooth), while on the other hand, for small pairing, we recover the sinx behaviour: 



(3.14) 



ev 



Ah,Ai,<<1 ttI 



—-ArAl sinx- 



(3.15) 



It is instructive to plot ( p. 13 ) in various regimes. Taking symmetric pairing Ar — Al = As to start with, we can 
see how Ishii's sawtooth is rounded off progressively to a sinx function as As is taken from 1 (perfect Andreev) to 
smaller and smaller values (i.e. for progressively more normal reflection at the contacts). This is illustrated in figure 




0.5 



1 . 5 



2 . 5 



3 



FIG. 5. Josephson current in the symmetric junction as a function of x, in units of e/ttI. The five plots are respectively for 
\AB\/t = 1,0.9,0.8,0.6 and 0.3. 



It is important to note that the expression ( 3.13 ) for the Josephson current is valid for independent arbitrary values 
of the boundary pairings, and thus covers the case of asymmetric junctions already studied for example in [ |l^. In 
this work, the shape of the current-phase relationship was still the sawtooth function, with critical current depending 
on the asymmetry between the pairings. The sawtooth result implies that effective perfect Andreev conditions were 
imposed, and thus that no normal reflection occurred at the conta cts. O ur expression for the current in a long junction 
thus covers a wider regime than the one in [ When plotting ( [3.13 ) for various asymmetries, the graphs look very 
similar to flgure |[ 

The formidable looking expression (3.13) can be considerably simplifled when A^ ~ Ar « 1. We first state the 



approximations and then justify them afterwards. We can approximate: 



arccos 



5(cosx) 



(.g2(l) + A2«|,)i/2 



(3.16) 



We may Taylor expand the x dependence of the other factors near |xl ~ since they give essentially a constant 
elsewhere. Thus: 



g(cos2 I) «4(1-A^)(1-A0 



(3.17) 
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giving: 



IXI < 



(3.18) 



Note that the last factor vanishes at |x| ^ but approaches 1 for tt — |x| >> |2 — A/j — Al|. For small |2 — A/j — A^j 
we find that the maximum current occurs at: 

2/3 

TT- \Xm\ « 7^1/3 



2(2- Ai- Afl) 



and has a value: 



1 



2(2- 



-1 2/3' 



Now let us consider the justification for these approximations. First note that: 
where, for convenience, we have defined: 

e = (2 - Ai - An)/vF. 

Thus, 

g(cosx) 
_(g2(i) + A24)i/2 

except near |x| « tt where arccos becomes singular, behaving as: 

arccos(— 1 + 6) ~ tt — V26 

Thus, near \x\ « tt we may write: 



arccos 



g{cosx) 



(g2(l) + A24)i/2 



n~y/{n-\x\r + 0{e^). 



Noting that, at the maximum xm, (tt — |x|) oc e^^^ we see that in this range of 



X- 



arccos 



Kcosx) 



(g2(i) + A24)i/2 



x[l + Oie'/% 



(3.19) 
(3.20) 

(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 

(3.26) 



Thus our simple a ppro ximation to arccos is everywhere valid. The other corrections from expanding sinx in the 
numerator in Eq. ( |3.13 ) and the expression in the denominator give multiplicative corrections of 0(e) or 0{tt — |x|)^ 
which is ©(e'*/^) near xm- 

In the case of perfect Andreev reflection at both boundaries, we can in fact solve for the ground state energy (and 
thus the Josephson current) exactly for arbitrary junction length I and filling. Putting An = Al = t{= 1) in (3.3) 
directly yields after simple algebra the allowed quasiparticle momenta for the two sets of Andreev levels. We find 

(3.27) 



l-l 



The ground state energy is then again the sum of the energies of all occupied levels, which now becomes a simple 
geometric progression: 

N sin COS "'^^+^^ 

i?o = -2f^[cosa++cosa_] = -4t Hl!IlL_^(fzlL cos — ^ . (|xl < tt) (3.28) 

The Josephson current is then, for arbitrary length / and occupation number N ^ 



Hx) 



ttN -it{N+1) 

Aet 2(ili) cos 2(,_i) 



X 



I - 1 



2{i-iy 



\x\ < V 



(3.29) 



One can explicitly check that the large / limit reproduces (3.14). 
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IV. RENORMALIZATION GROUP ANALYSIS OF THE INTERACTING CASE 



In a standard way the bulk Hamiltonian can be approximated in the continuum hmit, valid at low energies, by a 
quantum field theory, corresponding to a Tomonaga-Luttinger liquid describing gapless charge and spin bosons. The 
first step is to write the lattice fermion operators in terms of left and right moving continuum fermion operators: 



-ikpx 



(4.1) 



where kp is the Fermi wave-vector and V'l,_r are assumed to vary slowly on the scale of the lattice spacing (which 
is set to 1). The resulting continuum Hamiltonian is then bosonized in terms of charge and spin boson, (pes with 
associated velocities, Vc.s and compactification radii Rc,s- We will normally set these velocity parameters to 1. The 
continuum fermion fields are written: 



exp 



exp 



— - — \- nRcdc ± „ ± TrRscbs 
2Rc 2Rs 



— ~ nRc<j)c ± ^ T TrRs4>s 



(4.2) 



(4.3) 



Here the bosons and dual bosons are written in terms of left and right moving components as: 

0(t, x) ^ 4>L{t + x) + 4>ii(t - x), 4^{t, x) ^ (j)L{t + x) ~ (t)R{t - x). 

The Lagrangian has conventional normalization: 

£ - (1/2) [d^cj^cd^c^c] + (1/2) [d^(Psd^(l>s\ ■ (4.4) 

Here we follow the conventions of [ |l^. Unfortunately, various other bosonization conventions are frequently used. 
In particular, the compactification radii, Res, which depend on the bulk interactions, are often removed from the 
bosonization formulae by rescaling the bosons, resulting in an unconventional normalization of the two terms in the 
Lagrangian. The resulting normalization constants are sometimes called gp^a- The relationship between parameters 
gp^a- used in [m|,^,0, the parameters Kp^^^ used in [||] and our parameters is: 



1 



2K, 



(4.5) 



p,<j 



In the case of SU{2) symmetry, Rg = For replusive bulk interactions R^ > l/\/27r and for attractive bulk 



interactions Rc < l/v27r. In the case of attractive interactions, there may be a gap for spin excitations depending on 
the detailed form of the bulk interactions. This occurs, for example, for the attractive Hubbard model. The presence 
of a bulk spin gap makes very little difference to our analysis. Essentially, we may just drop 0s from our formulas. 
We consider both cases below. P-, 
Free boundary conditions, which occur for = 0, correspond to:Ej 



^i(0) = e-'Vfl(0), 



(4.6) 



where the phase 9 depends on the boundary scattering potential Vb ■ Note that these boundary conditions correspond 
to only normal reflection, and thus zero Andreev reflection, so in this context it is appropriate to refer to them as 
"normal" b.c.'s. In bosonized form these b.c.'s become: 



(4.7) 



It is crucial to realise that these equations imply that (j)Rc,s may be regarded as the analytic continuation of (/)z,c,s to 
the negative x-axis: 



In particular, this implies 



4>cr{x) = -~4>cl{-x) + Rc9, (t>sR{x) = -4>sl{-x), {x > 0). 



(t>c{Q)~^24>Lc{0)- RcO. 



(4.8) 



(4.9) 



We now wish to consider the bosonized form of the boundary scattering potential, cx Vb and boundary pairing 
interaction, cx As in Eq. (1.7). The scattering potential is proportional to dxfj^ci^)- This has scaling dimension 1 and 
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hence is marginal. Note that boundary interactions are relevant if their dimension is d < 1 and irrelevant if d > 1. 
This is different than for bulk interactions in (1+1) dimensions due to the fact that boundary interactions are only 
integrated over time, not space. Thus a boundary scattering potential leads to a line of fixed points, characterized by 
a phase shift. 

On the other hand, the boundary pairing interaction leads to the term: 

Ase"''^LaV'fl/3 + h.c. cx |Ab| sin[2^i?e^,(0) + x] cos[(/.,(0)/i?,] ^ As cos[47ri?,(/.i,(0) - 2ttR19 + x], 

(4.10) 



where the b.c. of Eq. (4.7) was used in the last step and x is the phase of A^. This operator has dimension 2ttR^ 
and hence is irrelevant for repulsive bulk interactions but relevant for attractive bulk interactions. Thus we reach 
the important conclusion that a weak boundary pairing interaction. As, becomes progressively less important as 
T — !■ in the case of repulsive bulk interactions. This implies that the effective coupling of the superconductor to the 
Luttinger liquid, t" , renormalizes to since A^ oc {t")'^. 

In the case of attractive bulk interactions, Rc < 1/^/2tt the free boundary condition is an unstable fixed point. 
The "obvious" guess is that the Hamiltonian renormalizes to a boundary fixed point corresponding to the boundary 
condition 

^c(0) = -(x + 7r/2)/27ri?„ (/..(O) = 0. (4.11) 

Note that the boundary condition on the spin boson is unchanged. This is surely a reasonable assumption since the 
boundary interaction doesn't involve the spin boson. In fact this boundary condition is fixed by SU(2) symmetry. 
On the other hand, we are assuming that the effect of the relevant boundary sine-Gordon interaction is to pin the 
dual charge boson, (j)c{0), corresponding to a semi-classical anabjsis-cif the interaction at large A^. We note that the 
analogous assumption has been made in several other contexts.Boll3 It is generally believed that only Dirichlet and 
Neumann fixed points occur in the boundary sine-Gordon model for generic compactification radius. [We note that 
(j) = constant corresponds to a Dirichlet b.c. and = constant to a Neumann b.c. using the fact that d4>/dt — d(t>/dx] 
In order to shed more insight on this assumption, we discuss, in the next section, a different boundary model which 
is equivalent to this one under an exact duality transformation. In the present context this boundary condition 
corresponds to perfect Andreev reflection since it follows from Eq. ( p. 12 ). 



The consistency of this assumption can be checked by considering the renormalization group stability of the Andreev 
b.c. Note that the scaling dimension of boundary operators are different at this fixed point where we must use: 

^cr{x) = <t>cL{-x) + (x + 7r/2)/27ri?„ 0,fl,(x) = ~4,sl{-x), {x > 0). (4.12) 

In this case a further boundary pairing interaction is marginal, corresponding to shifting the condensate phase, x- 
The potentially relevant interaction corresponds to normal scattering. This corresponds to adding a term 

6H = VN^l„iO)i'B.am + h.c. (4.13) 

to the Hamiltonian obeying the Andreev b.c. Using the Andreev boundary condition, and letting 9 be the phase of 
Vn, this term reduces to: 

SH oc -\Vn\ sin[(/.,/i?, + e]oc -\VN\sm[2cbcL/ Rc + (x + 7r/2)/2TTRl + 6], (4.14) 

of dimension l/(27ri?^). This is irrelevant for attractive bulk interactions, but relevant in the repulsive case. Thus 
we conclude that the Andreev b.c. represents an attractive fixed point with attractive bulk interactions, so that our 
assumption that a boundary pairing interaction leads to a flow to the Andreev b.c. for Rc < 1/V2tt is consistent. 
On the other hand, in the case of repulsive bulk interactions we expect an RG flow from the Andreev fixed point to 
the normal fixed point. For the non-interacting case, both fixed points are marginal and no renormalization occurs. 
There is a line of fixed points connecting normal and Andreev fixed points along which the ratio of normal to Andreev 
scattering variesxontinously. This behaviour is very analogous to the backscattering problem for a single impurity in 
a quantum wire.Q 

The behaviour of the Josephson current in the case of attractive bulk interactions is especially interesting. From 
Sec. 3 we see, that for the non-interacting case, the Josephson current is a sawtooth function of amplitude evp/l 
when the boundary terms are fine-tuned to give perfect Andreev scattering at the Fermi surface. Otherwise, /(x) is a 
smooth function. As shown by Maslov et aLd, /(x) is also a sawtooth function in the presence of bulk interactions if 
the Andreev b.c. is applied to the bosonized theory, with the amplitude replaced by evp / {2ttR^I). The sawtooth form 
is a universal property of the Andreev fixed point. This universality of the 0(1/0 term in the groundstate energy is a 
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familiar aspect of conformal field theory. In cases where the boundary parameters are not fine-tuned to the Andreev 
boundary condition, but instead the Hamiltonian renormalizes to the Andreev fixed point we expect to recover the 
same sawtooth form of J(x) in the limit I oo. However, the finite size corrections will smooth out I{x) since a finite 
I cuts off the RG flow of the junction parameter A^. We expect that for large enough I we may use the expression 



Eq. (3.18) with 1 — A replaced by an effective value at scale I. We may identify: 

1 - As « Vn, (4.15) 
where Vn is the effe ctive normal scattering interaction introduced in our discussion of the RG stability of the Andreev 



fixed point in Eq. (4.13). This identification is reasonable since it can be checked, for the non- interacting case, that 



the normal reflection amplitude vanishes linearly in 1 — |Ab|. Thus we expect that: 



« (4-16) 



Substituting this expression into Eq. ( 3.18D gives an approximate expression for the Josephson current. As / increases, 



/(x) becomes a more and more rapidly varying function near |x| « tt. The maximum of J(x) occurs at: 

IXmI ^2,l/(2^^g)_l]/3 (4.17) 



and the critical current scales with / as: 



2TTm 



constant 



p[l/(27ri??)-l]/3 



(4.18) 



The smoothing of /(x) for finite size junctions due to renormalization effects is quite distinct from the finite size 
effects that occur in the non-interacting case, discussed in Sec. IIL These are suppressed by powers of 1/P and do 
not smooth out the sawtooth structure of /(x)- In particular the maximum remains at |xm| = 

We note that our result for the behaviour of the Josephson current near the Andreev fixed point is very different 
from that obtained in [ ^ although both treatments use bosonization and RG arguments. The difference arises in 
part because we take into account the Andreev b.c. in calculating the RG scaling of the normal reflection amplitude 
and in part because we take into account the singular dependence of the Josephson current on the normal reflection 
amplitude. 

In the case of repulsive bulk interactions and almost perfectly fine-tuned junction parameters a flow away from the 
Andreev b.c. occurs with increasing junction length. In this case the effective parameter 1 — |Ab(Z)| increases with 
increasing I so that the sawtooth singularity is smoothed out as I increases. 

As mentioned above, in the case of attractive bulk interactions a spin gap sometimes occurs, for example in the 
[/ < Hubbard model. This has essentially no effect on the boundary RG discussed above since the spin boson didn't 
play any role. Essentially the spin boson is assumed to always obey the Dirichelt b.c, 0s (0) — throughout the RG 
flow which only affects the b.c.'s on the charge boson. In the case where there is a spin gap, (l)s{x) is pinned at all 
points in space; this is completely compatible with the assumption about the b.c. 

We find the flow to the Andreev b.c. in the attractive case especially remarkable because, as explained in the previous 
section, in the non-interacting case perfect Andreev scattering can only be achieved by flne-tuning parameters. Thus, 
in the interacting case, the RG flow must "find" the special value of the parameters at which the normal scattering 
vanishes. 



V. DUALITY TRANSFORMATION 



In an effort to make more plausible the conjectures about RG fiows in the previous section and in order to make 
contact with previous work on quantum impurity problems we present in this section an exact duality transformation 
from the lattice boundary pairing model with bulk interactions of the previous-jSection to a lattice model with both 
bulk and boundary magnetic fields. This is related to the previously studiedcJ S=l/2 xxz chain with a transverse 
boundary field. These latter models are perhaps easier to understand intuitively because semi-classical approximations 
hold to some extent. Furthermore there is an instructive difference between the Hubbard chain and pure spin chain 
corresponding to a sort of breakdown of "spin-charge separati on" for strong boundary fields. 

We begin with the (semi- infinite) boundary pairing of Eq. (1^) with the Hubbard interaction of Eq. (|l.9| ) added. 
We then apply the well-known duality transformation which changes the sign of the Hubbard coupling constant, J7, 
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and interchanges charge and spin operators. This is essentially a particle-hole transformation for spin up electrons 
only: 

^ (-irV'lT. i^jl^^jl- (5-1) 

This maps the hopping term into itself and the Hubbard interaction into (— l)x itself. The chemical potential term 
is mapped into a magnetic field in the z-direction: 

^l^ja^^-i^lC'^Ui^jP- (5.2) 

Thus a non-zero chemical potential, corresponding to average particle number < rij 1 maps into a non-zero bulk 
magnetic field in the z-direction. Note however, that the dual model has zero chemical potential so it remains at 
half-filling. Longer range density-density interactions map into z-z magnetic exchange interactions. The boundary 
scattering term, Vb maps into a modified boundary field in the z-direction. The boundary pairing interaction is 
mapped into: 

Hb^ A*s^pl^ijU+h.c. (5.3) 

This corresponds to a boundary magnetic field lying in the xy plane, transverse to the bulk field, of magnitude 2|Ab| 
and direction determined by the phase of . 

This dual model is especially easy to analyse in the case where C/ < so that there is a spin gap in the Hubbard 
model. The dual model, with U > and half-filling has a gap for charge excitations. The remaining gapless spin 
excitations are approximately described by the Heisenberg model with the appropriate magnetic fields. To make 
this correspondance more precise, when \U\ >> t, the correspondance holds for the lattice models with an effective 
Heisenberg exchange interaction /U . For smaller U the correspondance still holds for the low energy degrees of 
freedom. Even in situations where the original spin excitations were not gapped so that the dual charge excitations 
are not gapped, we might expect some sort of correspondance with the Heisenberg model at low energies due to 
spin-charge separation. 

The xxz S=l/2 spin model with a transverse boundary field (but no bulk field) was analysed in [ |l^. There it was 
shown that the bosonized version is the boundary sine-Gordon model with a boundary interaction which is relevant 
along the entire bulk xxz critical line and it was conjectured that an RG flow to the Neumann b.c. occurs. In the 
particular case of the xx model this can be proven exactly using Ising model duality transformations .Ej The semi- 
classical interpretation of the Neumann b.c. in this case is that the boundary spin is polarized in the direction of the 
boundary field. This analysis can be easily extended to include a bulk magnetic field in the z-direction. As shown in 
[ , the dimension of the transverse boundary field is 2'kB? where R is the compactification radius of the boson in 
the spin chain. [To fix our conventions, the transverse staggered correlation exponent is also 2t:R^] This becomes 
marginal for the isotropic xxx model but is relevant along the entire (zero field) xxz critical line. The depcndpice of 
the radius on a magnetic field applied to the Heisenberg model has been calculated from the Bethe ansatz.EZlE3 The 
effect is again to decrease the radius, hence making the transverse boundary field relevant. Thus it is again natural to 
conjecture a flow to a flxed boundary spin polarized in the direction of the boundary field. Thus, flow to the Andreev 
b.c. in the boundary pairing model is dual to flow to a fixed spin b.c. in the boundary field model. 

This analysis can be extended to more general models with longer range bulk interactions. In particular, we may 
consider cases in which the spin excitations are not gapped in the original model so that charge excitations are not 
gapped in the dual model. Again it seems plausible that even a weak transverse boundary field produces a fiow to a 
polarized spin boundary fixed point. However, we now encounter another interesting phenomenon. If the boundary 
field is too strong it suppresses this RG flow. This can be seen from the fact that, in the limit of a very strong 
transverse boundary field, one electron gets trapped on the first site with probability 1 in a state with spin polarized 
along the transverse field direction. Since the hopping term adds or removes an electron from site 1 it produces a 
high energy state, with energy of order the boundary field, |Ab|. All such processes are suppressed for \Ab\ » t 
meaning that the first site decouples from all the others which therefore obey a free b.c. Thus, in the finite U model 
the spin-polarized Neumann b.c. should not be thought of as occuring at infinite boundary field, but rather at a finite 
value. On the other hand, in the Heisenberg model we may indeed think of the spin polarized fixed point as occurring 
at infinite boundary field since the magnetic exchange interaction isn't suppressed by the strong field. Thus we see 
that the limit U ^ oo and |Ab| ^ oo do not commute. 

The above observation provides another way of understanding the perhaps surprising discovery in the previous 
sections that the Andreev fixed point does not occur at oo boundary pairing strength but rather at a fine-tuned finite 
value. In this model at very strong Ab we may think of a sort of Andreev boundstate occurring on the first site 
corresponding to a linear combination of the vacuum and filled state: 

|0>+e*>^|t,i>. (5.4) 
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Since the hopping term always turns this Andreev boundstate into a state with a single electron at site 1 it produces a 
high energy state and its effects are therefore suppressed when \Ab\ >> t. In the original SN interface model we may 
think of the Andreev boundstate as blocking electron transport across the interface and hence suppressing Andreev 
scattering. 

This research was supported by NSERC of Canada and the NSF under grant No. PHY94-07194. 



APPENDIX A: INTEGRATING OUT THE SUPERCONDUCTING ELECTRONS 

In this appendix, we outline the steps leading to the correspondence between the parameters in the interface and 
boundary models. 

Let us thus consider the interface model, which has nonzero gap on sites j < 0. Our strategy will consist in 
integrating these out. For simplicity, let us use the notation tp^^ for the fields on site one. Omitting all sites with 
j > 1 for ease of notation, we can write down the contribution to the imaginary-time action coming from the gapped 
side and its coupling to the boundary fields: 

H'Vf,\u;)^o.M + /i.c.)}. (Al) 
Fourier transforming as Vja = f Jq dksmk{j — l)ip„{k) (for j < 0) and using the Bogoliubov transformation 

(A2) 



V^t(w,/c) \ _ f u{k) -v*{k)\f ri+{uj,k) 

i^\{-u^,k) ) - \v{k) u*{k) J{nU-^,k) 



where 



and e{k) = ~2t' cos k — n and E{k) = ^e^(fc) + |Ap, we arrive at the form 

S^ = ^Y.l I dk [rjliLO,k)[iiJ + E{k)]r]^{iJ,k)[-t''r]l{u;,k)[u*{k)i;^^^^^ 

-t"r]l {u, k) [u* (fc) Vf^ (w) + V* (fc) ^ (-w)] sin fc + /i.e.] } . (A4) 
Integrating out the r] fields finally gives the boundary action 

^B = \ E {[^^^i^'^) - C2{u:)]i,^J + (ci(a;)A^f/(u;)Vf/(-u;) + h.c.)] (A5) 

where the coefficients Cj(a;), appearing respectively in front of the pairing-like and potential-like amplitudes, are given 
by 

-\/a;^-|-|A|^ -I- |A|^ 

in which 



M = s^-2t' + /z + iVw2 + |^|2 + ^2t' + /X + i + |A|2. (A7) 
The desired correspondence between the interface and boundary parameters thus takes the form (when a; ^ 0) 

As = Aci(a;^0), Vb = F - 02(0; ^ 0). (A8) 
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Taking the limit explicitly reproduces equations (2.18). The frequency-dependent terms are suppressed by powers of 
w/jAj and have thus been ignored at energies well below the gap. Furthermore, in the presence of interactions, we 
expect the above procedure to work as well, namely that the final result is simply some effective boundary pairing 
and scattering potentials. 
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